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^This  report  describes  the  new  portable  platform  and  gravity  meter 
system  which  has  been  assembled  at  the  Woods  Hole  Oceanographic  Institution. 
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The  new  platform  represents  a  major  reduction  in  both  size  and  weight  over 
other  platforms  suitable  for  gravity  measurement.  The  second  major  part  of 
this  system  is  a  new  gravity  readout  which  interfaces  with  the  VSA,  process¬ 
es  the  VSA  output,  and  prepares  the  resulting  filtered  acceleration  data 
for  output  to  the  acquisition  system.  The  readout  has  been  designed  to  allow 
flexible  use  of  the  gravity  system  on  a  variety  of  vehicles,  including  ships, 
submarines  and  aircraft.  The  third  part  of  this  new  meter  is  the  data  ac¬ 
quisition  system.  It  consists  of  a  microprocessor  interfaced  to  a  Kennedy 
9-track  tape  drive.  Both  the  platform  and  the  readout  are  connected  to  the 
microprocessor. 

Results  are  presented  from  Endeavor  cruise  88  that  demonstrate  the 
ability  of  the  platform  to  stabilize  the  gravity  meter  and  for  the  gravity 
system  to  produce  raw  data  with  a  resolution  of  48  milligals  at  a  sampling 
rate  of  10  Hz.  44Kgital  signal  processing  techniques  which  were  used  to  fil¬ 
ter  the  data  and  extract  the  gravity  signal  with  a  resolution  of  0.48  milli¬ 
gals  are  also  discussed. 
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ABSTRACT 

This  report  describes  the  new  portable  platform  and  gravity  meter  system 
trtilch  has  been  assembled  at  the  Woods  Hole  Oceanographic  Institution.  It 
consists  of  three  functionally  distinct  parts. 

The  first  of  these  Is  a  recently  developed  gyro -stab 111 zed  two-axis  plat¬ 
form.  This  platform  has  been  designed  to  carry  the  vibrating  string  accelero¬ 
meter  (VSA)  and  Its  associated  oven  assembly  as  the  gravity  sensor.  The  new 
platform  represents  a  major  reduction  In  both  size  and  weight  over  other 
platforms  suitable  for  gravity  measurement.  The  second  major  part  of  this 
system  Is  a  new  gravity  readout  which  Interfaces  with  the  VSA,  processes  the 
VSA  output,  and  prepares  the  resulting  filtered  acceleration  data  for  output 
to  the  acquisition  system.  The  readout  has  been  designed  to  allow  flexible 
use  of  the  gravity  system  on  a  variety  of  vehicles,  Including  ships,  submar¬ 
ines  and  aircraft.  The  third  part  of  this  new  meter  Is  the  data  acquisition 
system.  It  consists  of  a  microprocessor  Interfaced  to  a  Kennedy  9-track  tape 
drive.  Both  the  platform  and  the  readout  are  connected  to  the  microprocessor . 

Results  are  presented  from  Endeavor  cruise  88  that  demonstrate  the  ability 
of  the  platform  to  stabilize  the  gravity  meter  and  for  the  gravity  system  to 
produce  raw  data  with  a  resolution  of  48  mllllgals  at  a  sampling  rate  of  10 
Hz.  Digital  signal  processing  techniques  which  were  used  to  filter  the  data 
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and  extract  the  gravity  signal  with  a  resolution  of  0.48  mllllgals  are  also 
discussed. 

SYSTEM  DESIGN  OBJECTIVES 

Primary  design  objective  for  tile  new  gravity  platform  Is  to  provide  a 
small  two-axis  stable  table  which  can  be  Installed  In  a  variety  of  vehicles 
and  which  will  carry  the  existing  VSA  gravity  sensor.  The  range  of  vehicles 
that  Is  considered  acceptable  for  gravity  work  Includes  surface  ships,  submar¬ 
ines  and  airplanes.  To  take  advantage  of  the  opportunity  provided  by  such  a 
selection.  It  Is  necessary  that  the  system  be  easily  portable.  Since  some  of 
these  vehicles  will  have  little  space  to  devote  to  additional  Instrumentation, 
the  platform  shall  be  as  small  and  compact  as  possible. 

The  platform  should  maintain  level  with  sufficient  accuracy  that  gravity 
data  can  be  resolved  to  1  mill 1 gal  while  the  system  is  In  operation  on  a  sur- 


face  ship  In  moderate  seas.  The  stable  platform  Inertial  data,  which  consists 
of  the  Instantaneous  accelerometer  outputs  and  the  Integrated  accelerations, 
or  Inertial  velocities,  are  provided  as  outputs  for  subsequent  analysis  of  the 
vehicle's  dynamic  motion.  Glmbal  angle  resolvers  provide  output  of  pitch  and 
roll  data. 

Data  Input  to  the  platform  are  speed,  heading  and  latitude.  Because  It  Is 
anticipated  that  the  stable  platform  will  be  used  In  airplanes.  It  must  be 
capable  of  handling  speeds  of  several  hundred  knots  as  well  as  the  more 
leisurely  pace  of  the  surface  ship  or  submarine. 

We  were  fortunate  to  have  received  from  the  Charles  Stark  Draper  Labora¬ 
tories  { CSDL )  two  engineering  models  of  the  Apollo  Inertial  navigation  plat¬ 
form.  These  were  given  to  us  as  a  source  of  gyros  and  accelerometers  and 
other  parts  that  could  be  fit  Into  our  design.  From  these  Apollo  navigators 
we  obtained  two  gyros,  two  accelerometers  and  four  torque  motorbearing 
assemblies  with  sllprlngs.  Using  these  parts  and  the  design  goals  stated 
above,  we  developed  the  design  for  the  small  stable  table. 

Platform  Description 

The  platform  that  was  designed  to  meet  these  requirements  Is  depicted  In 
simplified  schematic  form  in  figure  1.  This  figure  shows  only  one  axis  of  the 
two,  since  they  are  identical.  The  electronics  and  the  computer  are  mounted 
separately  from  the  gimbals,  to  which  they  are  connected  by  a  10  ft.  cable. 

Because  the  Apollo  accelerometers  are  pulse-torqued,  it  was  decided  that 
the  best  approach  for  the  electronic  design  would  be  to  use  a  digital  control 
loop.  Included  within  the  loop  would  be  a  small  computer.  In  this  way,  data 
from  accelerometers  and  navigation  devices  would  be  used  to  calculate  directly 
the  torqulng  rates  applied  to  the  gyros.  Software  would  replace  the  mechanical 
or  electrical  Integrators  that  have  been  traditionally  used,  and  the  system 
could  be  easily  adapted  to  different  configurations. 

Consultation  for  this  design  was  provided  by  A.  Gulovson  from  CSDL.  With 
his  help,  we  developed  the  platform  operating  requirements  and  the  digital 
control  loop  concept.  Based  on  his  past  experiences  with  Inertial  navigation 
systems  and  his  associations  with  the  people  who  were  Involved  with  the  Apollo 
system,  he  was  able  to  provide  us  with  circuit  diagrams  and  operating 
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procedures  for  the  major  portions  of  the  system.  These  Include: 

1)  Gyro  and  accelerometer  torqulng  loop 

2)  Demodulators  and  servo  amplifiers 

3)  Resolver  demodulators 

4)  Modlng  logic 

5)  Heater  control  circuits 

6)  Oscillator  and  timing  circuits 

7)  AC  power  supplies 

In  addition  to  the  circuit  diagrams  and  the  Apollo  specification  drawings, 
Gulovson  provided  us  with  a  computer  program  which  he  used  to  test  the  concept 
of  the  digital  control  loop.  Within  this  program  were  the  algorithms  and 
equations  which  we  Included  In  our  real-time  computer  program. 

MAJOR  SYSTEM  COMPONENTS  AND  OPERATION 

The  central  elements  of  an  inertial  platform  are  the  accelerometers  and 
gyros.  These  components  determine  the  maximum  performance  which  can  be 
attained  from  the  system.  The  gyro  Is  used  for  platform  position  control, 
while  the  accelerometer  Is  used  for  vertical  reference.  These  Instruments  are 
mounted  on  the  stable  member,  which  is  in  turn  supported  by  a  gimbaled  system. 

The  short-term  vertical  reference  Is  provided  by  the  Inertial  Rate  Inte¬ 
grating  Gyro.  The  unit  consists  of  a  float,  or  housing,  for  the  gyro  motor,  a 
signal  generator  which  indicates  float  position  relative  to  a  null  point,  and 
a  torqulng  mechanism  which  can  rotate  the  float  about  the  input  axis.  The 
float  Is  heated  to  135  degrees  F.  and  maintained  at  this  temperature  so  that 
It  Is  neutrally  buoyant  within  a  fluid  filled  case.  It  Is  loosely  held  near 
the  center  of  the  case  by  an  electro -magnetic  suspension  system.  A  two-  phase 
motor  at  the  center  of  the  float  turns  the  gyro  wheel  at  24000  rpm.  Position 
of  the  float  about  a  null  point  Is  sensed  by  the  signal  generator  which  has 
output  of  11.8  mllllvolts/mlcro-radlan.  A  torqulng  command  signal  may  be  Input 
to  the  gyro  to  rotate  It  about  Its  Input  axis  in  order  to  compensate  for  such 
things  as  gyro  drift,  earth  rate  and  velocity. 

The  long-term  vertical  reference  Is  provided  by  the  Pulse  Integrating 
Pendulous  Accelerometer  (PIP).  This  has  a  heated  float  element  similar  to 
that  of  the  gyro.  The  float  Is  weighted  and  suspended,  creating  a  pendulum 
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which  rotates  freely  about  one  axis.  When  the  PIP  is  accelerated,  the  float 
rotates  by  an  amount  proportional  to  the  acceleration.  This  rotation  away 
from  the  null  point  is  sensed  by  the  signal  generator,  with  the  phase  of  the 
ilgnal  Indicating  on  which  side  of  the  null  the  pendulum  has  moved.  The  phase 
Information  from  the  signal  generator  Is  amplified  and  output  to  a  torque 
generator  attached  to  the  float.  Utilizing  this  data,  the  torque  generator 
returns  the  float  to  the  null  point.  Thus,  as  continuous  pulses  of  current 
flow  Into  the  torque  generator,  the  float  Is  kept  hovering  about  zero.  These 
pulses,  when  summed,  represent  the  amount  of  restoring  force  required  to 
maintain  the  float  at  null.  This  Is  directly  proportional  to  the  acceleration 
along  the  sensitive  axis  of  the  Instrument. 

Glmbal  Rings 

The  stable  member  of  the  gravity  platform  Is  held  within  the  framework  of 
a  two-axis  glmbal  system.  This  stable  member  Is  made  from  a  0.5  Inch  hardened 
aluminium  plate  approximately  12  X  14  Inches,  with  a  cutout  to  allow  the  VSA 
oven  to  be  suspended  through  It.  The  gyros  are  mounted  In  a  two-piece  fixture 
that  places  one  above  and  one  below  the  aluminium  plate.  The  accelerometers 
are  housed  In  a  one-piece  block  and  also  mounted  above  and  below  the  plate. 
Symmetry  and  distribution  of  mass  In  tile  platform  design  have  been  maintained 
to  keep  the  center  of  mass  as  close  as  possible  to  the  axis  center. 
Preamplifier  circuit  cards  for  the  temperature  sensors  and  the  signal 
generators  are  mounted  In  edge  connectors  near  each  Instrument. 

The  stable  member  Is  supported  by  a  pair  of  torque  motor  assemblies.  The 
weight  Is  carried  by  stub  shafts  riding  In  a  set  of  bearings.  These  bearings 
transfer  the  load  to  the  outer  case  of  the  torque  motor.  A  face  plate  Is 
fastened  to  the  Inner  end  of  the  stub  shafts  and  this  provides  the  mounting 
surface  for  the  stable  member.  The  outer  cases  of  the  torque  motors  are 
mounted  In  fixtures  that  are  Integral  components  of  the  glmbal  rings.  Thus, 
the  inner  ring  consists  of  two  loops  of  formed  aluminium  tubing  which 
Interconnects  the  fixtures  for  holding  the  torque -mo tors.  An  outer  ring  made 
from  three  loops  of  tubing  plus  the  torque-motor  housings  provides  support  for 
the  Inner  ring.  The  outer  ring  Is  shock-mounted  and  attached  to  a  baseplate 
and  support  structure.  Pitch-axis  control  Is  provided  by  rotation  of  the 
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stable  member  on  the  torque  motor  axis  Inside  the  Inner  ring.  The  roll-axis 
pivots  on  the  axis  of  the  torquers  fastened  to  the  outer  glmbal  ring.  Both 
axes  can  rotate  through  360  degrees. 

Besides  supplying  structural  support  between  the  two  glmbal  rings  and  the 
stable  member,  the  torquer  assemblies  provide  for  monitoring  of  the  glmbal 
angle  and  for  the  passage  of  electrical  signals  between  rings.  Each  unit 
also  contains  a  DC  torque  motor  used  to  maintain  the  position  of  the  glmbal 
controlled  when  an  external  force  acts  to  rotate  it  out  of  alignment.  Each 
pair  of  torqu^  -semblies  has  a  16X  and  a  IX  resolver  for  angle  readout  and 
manual  glmbal  control  in  the  "standby  mode."  The  IX  resolvers  are  connected 
to  a  12-bit  synchro-to-dlgltal  converter  so  that  the  pitch  and  roll  angles  are 
available  to  the  data  acquisition  system  with  a  resolution  of  8.5 
arc-minutes.  The  center  of  each  assembly  contains  a  slipring  unit  with  48 
Individual  circuits.  All  electrical  signals.  Including  shields  and  ground 
pass  to  the  stable  member  through  the  sllprlngs. 

Platform  Concepts 

There  are  two  major  control  loops  within  the  stable  platform  system.  The 
first  of  these  Is  for  glmbal  stabilization.  This  loop  connects  the  output 
from  the  signal  generator  of  each  gyro  to  demodulators  and  amplifiers  which 
drive  the  torque  motors  on  the  roll  and  pitch  platform  axis.  Since  the  gyro 
Is  mounted  rigidly  to  the  platform,  the  loop  Is  closed  through  mechanical 
coupling  between  them.  If  the  platform  Is  rotated  by  some  outside  force,  the 
motion  Is  sensed  by  the  gyro  as  a  misalignment  between  the  float  and  its  null 
point.  The  signal  generator  output  drives  the  torque  motors  restoring  the 
platform  to  Its  prior  position.  This  loop  has  sufficient  speed  to  maintain 
very  tight  coupling  between  the  gimbals  and  the  gyro.  Platform  errors  due  to 
the  stabilization  loop  are  less  than  15  arc-seconds. 

The  second  major  loop  Is  the  Vertical  Erection  Loop.  This  feedback  system 
senses  the  local  vertical,  and  torques  the  gyros  to  align  with  a  plane  which 
Is  perpendicular  to  It.  In  operation,  data  from  the  accel erometers  are 
processed  In  the  computer,  along  with  navigation  Information  so  that  a  level 
plane  Is  maintained  by  the  generation  of  gyro  torqulng  rates  which  compensate 
for  the  earth's  rotation  and  for  vehicle  translation  and  acceleration. 
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In  the  new  gravity  platform,  the  Vertical  Erection  Loop  Is  a  digital 
system.  A  computer  with  a  real  time  data  processing  program  (Fig.  2)  is 
continuously  being  updated  with  acceleration,  speed,  heading,  and  latitude. 
With  these  data,  gyro  torqulng  rates  are  calculated.  These  rate  numbers  are 
then  converted  to  the  pulse  trains  which  are  Input  directly  to  the  gyro 
torqulng  electronics.  The  accelerometers  are  sampled  every  5  milliseconds.  A 
new  torque  command  is  issued  to  each  gyro  at  the  same  time.  The  acceleration 
samples  are  summed  for  one  second  to  generate  the  average  acceleration  value 
for  the  one  second  Interval.  This  value  becomes  the  current  input  to  the 
navigation  subroutine  (NAV79).  Speed,  heading,  and  latitude  from  the 
vehicle's  sensors  are  updated  as  available.  NAV79  converts  the  accelerations 
to  velocities,  filters  these  velocities  in  a  lag/lead  network  and  adds  a 
correction  for  earth  rate.  The  results  of  these  calculations  are  scaled  to 
produce  the  correct  torquing  rate  for  each  gyro. 

Systems  can  be  designed  which  will  theoretically  Indicate  the  vertical 
despite  movements  of  the  system  in  any  fashion  about  the  earth's  surface. 
However,  if  an  error  in  the  vertical  indication,  or  an  Initial  offset,  Is 
introduced,  these  theoretical  systems  will  oscillate  Indefinitely  at  a  period 
of  close  to  84  minutes.  Schuler  first  demonstrated  this  fact  for  a  compound 
pendulum  In  1923.  The  same  Is  true  for  all  modern  stable  table  designs.  The 
84  minute  period  Is  determined  by  purely  physical  quantities:  G,  earth 
gravity,  and  R,  earth  radius.  Since  errors  can  never  be  eliminated  In  real 
systems.  In  order  to  prevent  continuous  oscillation,  Insensitivity  to  vehicle 
motion  must  be  sacrificed.  A  damping  filter  is  Included  In  most  designs  which 
effectively  increases  the  natural  frequency  of  the  Instrument  above  the 
Schuler  frequency.  Normally,  this  damping  is  small  so  that  the  system  remains 
fairly  insensitive  to  motion,  but  oscillations  die  out  with  time.  In  some 
situations,  the  natural  frequency  Is  Increased  to  the  point  that  the  period  of 
damped  oscillation  Is  on  the  order  of  a  few  minutes.  The  latter  type  of 
filter  Is  used  particularly  at  startup  time  when  It  Is  desired  to  have  a 
system  alllgned  with  the  local  vertical  In  a  short  amount  of  time.  If  a 
system  continues  to  operate  In  a  highly  damped  "startup  mode",  It  will  go  off 
level  when  a  change  In  velocity  occurs,  but  the  instrument  will  recover  In  a 
comparatively  short  amount  of  time.  When  the  primary  purpose  of  a  system  Is 
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Figure  2  The  Real  Time  Computer  Program.  The  background  program  provides  acceleration  data  and 
outputs  gyro  torquing  commands  every  5  ms.;  every  second  a  new  torquing  rate  is 
calculated.  The  foreground  program  monitors  the  operational  mode  of  the  platform, 
provides  navigation  information,  updates  the  displays  and  communicates  with  the  data 
link . 


the  measurement  of  the  vertical  gravity  component,  there  Is  a  basic  tradeoff 
between  the  desire  for  Insensitivity  to  vehicle  motion,  and  the  need  for  quick 
recovery  from  an  off-level  state.  In  the  present  <Vc'gn,  damping  Is 
controlled  by  a  digital  filter.  Its  characteristics  can  be  modified  by  the 
operator,  thereby  allowing  for  experimentation  with  the  tuning  and  damping 
applied  to  the  gyro  torquing  loops.  It  Is  anticipated  that  the  platform  can 
be  adjusted  for  best  performance,  depending  on  the  particular  vehicle. 

A  digital  computer  inside  the  loop  of  an  Inertial  system  Is  not  a  new 
idea.  The  usual  implementation  has  been  to  use  counting  circuits  and  adders 
to  emulate  the  functions  of  an  analog  computer.  Analog  computers  were 
necessary  to  perform  the  multiple  Integrations  and  summations  for  conversion 
of  accelerations  Into  distance  and  to  maintain  the  proper  gyro  torquing 
rates.  However,  the  microprocessor  offers  new  opportunities  for  performing 
these  tasks.  In  this  system,  the  Integrations  and  filtering  are  performed 
using  digital  signal  processing  techniques.  The  result  Is  that  good  control 
is  maintained  over  the  performance  of  the  system.  The  software  Is  easily 
modified  or  adapted  to  Interface  with  specific  navigation  sensors  or 
environments. 

The  Microcomputer  System 

It  was  desired  to  keep  the  on-board  microcomputer  system  as  simple  as 
possible.  In  order  to  avoid  major  computer  hardware  development  work.  A  board- 
level  system  was  selected  to  allow  flexibility  in  hardware  conflguatlon  yet 
provide  us  with  well-defined  operational  equipment.  It  was  also  necessary  to 
maintain  software  compatibility  with  a  microcomputer  to  be  used  for  software 
development. 

The  microcomputer  selected  for  use  in  the  platform  electronics  consists  of 
a  collection  of  Standard  (STD)  Bus  microprocessor  components  manufactured  by 
the  Mostek  Corperatlon.  This  set  Includes  a  Z-80  CPU  card,  32k  RAM  memory 
card,  3  parallel  I/O  cards,  1  serial  I/O  card,  and  one  card  with  a  rudimentary 
operating  system  In  ROM  plus  a  TTY  port.  The  I/O  boards  Interface  the 
microcomputer  system  to  the  gyro- torquing  electronics,  accelerometers,  and 
other  sensors  and  displays  which  make  up  the  platform  electronics.  The  TTY 
port  on  tiie  ROM  card  Is  connected  to  a  larger  computer  system. 
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It  Is  the  Interface  to  the  larger  computer  system  which  lends  flexibility 
and  programming  ease  to  the  platform  computer.  The  larger  system  Is  also  a 
Z-80  CPU,  running  under  a  CP/M  operating  system.  With  this  configuration, 
software  modules  can  be  developed  using  the  higher-level  languages  and 
programming  aids  available  for  the  CP/M  system.  These  can  be  linked  to 
modules  written  in  the  Z-80  assembly  language  and  the  complete  program  loaded 
to  form  a  memory  Image  of  the  desired  computer  program  for  the  platform. 

Then,  using  a  program  written  for  the  purpose,  plus  the  absolute  loader  In  the 
Mostek  ROM,  the  memory  Image  Is  down-loaded  Into  the  RAM  in  the  platform 
computer  where  It  Is  executed.  This  arrangement  has  been  very  helpful  with 
software  development,  allowing  for  quick  turnaround  In  the  debugging  process. 

The  computer-to-computer  link  allows  for  convenient  changes  to  the 
software.  It  is  through  this  mechanism  that  the  operational  response  of  the 
platform  can  be  reconfigured  to  test  new  Ideas  or  to  match  performance  to 
different  vehicles. 

Additionally,  it  was  discovered  that  this  link  was  the  solution  to 
platform  data-acquisitlon.  By  modifying  the  FORTRAN  library  drivers,  a  WRITE 
statement  was  made  to  transmit  data  up  the  link  to  the  CP/M  computer.  This 
bi-directional  link  provides  the  ability  to  record  and  display  In  real  time 
any  of  the  numerical  results  of  NAV79  or  the  data  input  to  It.  From  these 
exchanges  of  Information,  the  behavior  of  the  platform  can  be  analyzed,  with 
modifications  made  to  system  constants  and  damping  parameters  while  the 
platform  Is  operating. 

CRUISE  RESULTS 

The  first  sea  trials  for  the  new  stable  platform-gravity  meter  system  have 
been  completed  on  Endeavor,  cruise  88.  This  cruise  was  one  of  a  series  for 
the  Warm  Core  Rings  Experiment,  at  the  edge  of  the  Gulf  Stream.  We  worked  in 
cooperation  with  T.  Joyce  of  the  Physical  Oceanography  department,  who  provided 
us  with  space  on  the  Endeavor  and  access  to  the  Loran  navigation  and  gyro 
heading  data  which  is  collected  by  his  Acoustic  Profiler  of  Ocean  Current 
(APOC)  system.  These  data  were  available  to  us  In  real  time  on  the  SAIL  loop 
(Serial  ASCII  Interface  Loop),  which  interconnected  the  current  profiler  with 
the  various  sensors.  The  stable  platform  was  connected  to  the  SAIL  loop  as  a 
listening  device. 
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The  first  three  days  of  the  cruise  were  used  to  complete  an  Initial 
software  system  for  data  acquisition.  Concurrently,  the  stable  platform 
computer  was  being  Interfaced  to  the  ship's  doppler  speed  equipment.  On 
August  6,  the  platform  was  operated  for  the  first  time  at  sea. 

The  glmbal  servos  became  unstable  Immediately,  oscillating  about  the  null 
point.  This  condition  was  determined  to  be  caused  by  two  factors:  1)  the 
temperature  of  the  platform  was  slightly  higher  and  more  uniform  than  It  had 
been  In  the  laboratory,  since  the  gimbals  were  now  covered  by  a  plastic 
housing,  and  2)  vibrations  were  present.  The  elevated  temperature  raised  the 
frequency  response  of  the  gyros.  Increasing  the  gain  bandwidth  product.  The 
vibration  provided  a  mechanical  high-frequency  input.  The  Increased 
high-frequency  response  from  the  gyros  changed  the  gain-phase  relationship, 
lowering  the  phase  margin  of  this  tuned  feedback  system,  and  causing  an 
Instabllty  at  higher  frequencies.  The  noise  and  vibration  provided  Inputs  in 
the  frequency  range  where  the  loop  was  not  unstable.  This  loop  was  probably 
marginally  unstable  during  tests  In  the  lab,  although  this  was  not  observed. 
The  glmble  servos  were  made  to  function  adequately  by  removing  the  cover  and 
adjusting  the  overall  gain  of  the  servo  amplifiers.  This  was  only  a  temporary 
fix,  and  the  system  needs  to  be  retuned  by  changing  the  gain  and  phase  delay 
of  the  servo  amplifiers. 

It  was  observed  that  while  In  the  "navigate  mode"  there  was  Insufficient 
damping,  so  that  the  platform  oscillated  at  the  Schuler  period,  see  fig  3. 
However,  the  oscillations  did  die  out  over  time. 

In  the  "fast  erect  mode"  the  system  Is  much  more  responsive.  The  Inertial 
velocities  are  accurate  within  0.5  kts  of  the  rough  estimates  of  ship's 
speed.  In  this  mode,  the  Integral  of  the  acceleration  data  tracks  velocity 
data  Input  from  the  Doppler  speed  meter.  In  effect,  the  closed  loop  error 
signal  becomes  the  difference  between  Integrated  acceleration  and  velocity 
data.  However,  the  platform  remains  very  nearly  level  over  the  long  term, 
since  the  integral  of  the  acceleration  Is  finite,  and  the  deviation  from 
velocity  Is  small.  A  profile  of  "fast  erect"  response  is  shown  In  figure  4. 

It  can  be  seen  from  the  figure  that  there  Is  a  one  to  one  correspondence 
between  the  velocity  vectors  derived  from  the  Doppler  speed  meter  (resolved 
with  the  gyro  compass)  and  the  velocity  vectors  generated  by  Integrating  the 
data  from  the  two  accelerometers. 
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Inertial  Velocity  Quick  Erect  Mode.  It  is  seen  that  the  integrated  accelerations 
match  closely  with  the  speed  and  heading  data.  The  abrupt  shift  is  a  course  change. 
The  length  of  the  plot  is  20  minutes. 


An  error  In  estimation  of  the  bandwidth  of  the  ship's  heave  resulted  In  a 
phase-lock  loop  design  for  the  gravity  readout  having  Insufficient  dynamic 
range  to  remain  In  lock  for  normal  sea  states.  Therefore,  vertical 
acceleration  data  were  obtained  only  for  short  periods  of  time  during  calm 
seas.  Data  were  obtained  directly  from  the  vertical  sensor  using  a  frequency 
counter  In  order  to  correctly  estimate  the  range  for  the  redesign  of  the 
phase-locked  loop. 

The  gravity  meter  data  recorded  during  a  brief  calm  period  have  been 
reduced  using  numerical  analysis  techniques  to  be  discussed  below. 

Data  acquisition  programs  were  developed  for  recording  the  Inertial  data 
and  dynamic  motion  of  the  ship  as  sensed  by  the  gyro -stab 111 zed  platform. 

These  time-series  data  were  merged  with  the  frequency  counts  from  the  gravity 
readout  and  the  resultant  time  series  was  recorded  on  nine-track  magnetic 
tape.  Twenty  tapes  were  obtained.  The  data  Include  platform  behavior  In  the 
"fast  erect  mode"  and  In  the  "navigate  mode.”  They  also  Include  performance 
for  various  damping  parameters.  Recordings  of  LORAN,  speed  and  heading,  and 
vertical  acceleration  are  also  available  for  analysis.  Information  during  one 
XBT  star  was  recorded  for  subsequent  comparison  with  the  APOC  data. 

Processing  of  VSA  data  from  New  Platform 

Data  from  the  Vibrating  String  Accelerometer,  which  is  used  as  the  gravity 
sensor  In  the  system  (Bowln  et  al .1972 ) ,  Is  encoded  as  a  pair  of  narrow  band 
signals  which  are  brought  out  to  a  readout  unit  for  preliminary  processing. 

The  difference  In  frequency  between  this  pair  of  signals  Is  proportional  to 
the  acceleration  experienced  by  the  VSA  along  Its  Input  axis.  This  axis  Is 
aligned  to  the  local  vertical  by  the  platform.  The  difference  frequency  Is  on 

the  order  of  ^65  Hertz,  when  the  VSA  Is  static  with  Its  input  axis  along 
the  vertical.  A  rough  measure  of  Its  scale  factor  Is  then  6/65  or  about  15000 
mllllgals/Hz,  where  G  Is  the  normal,  average  vertical  gravity  component 
(980000  mllllgals). 

The  difference  frequency  Is  obtained  In  the  readout  by  applying  both  VSA 
signals  to  a  dual  balanced  mixer  which,  In  effect,  multiplies  the  two  Inputs 
and  generates  an  output  made  up  of  sum  and  difference  frequency  components. 

The  sum  component  Is  subsequently  filtered  out  with  an  analog  low  pass  filter 
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(LPF).  In  the  present  application.  In  which  typical  vertical  accelerations 
will  be  no  more  than  .3G,  the  signal  at  the  output  of  the  LPF  Is  a  narrow  band 
signal  centered  at  fd  Hz.  The  problem  of  obtaining  vertical  acceleration 
and  gravity  information  from  this  output  Is  thus  essentially  a  problem  In 
narrowband  FM  demodulation. 

The  information  encoded  In  the  modulating  frequency  Is  obtained  by 
counting  zero  crossings,  a  simple  technique  which  can  readily  be  used  for 
narrowband  FM  (Stark  and  Tuteur,  1979).  The  application  of  digital  signal 
processing  techniques  to  data  In  which  Information  Is  encoded  In  the  zero 
crossings  of  a  signal  Is  not  a  problem  often  encountered  In  the  related 
literature.  One  aspect  of  the  problem  Is  that  the  fundamental  sampling  rate  of 
the  Information  varies  with  the  timing  of  the  zero  crossings.  If  an  estimate 
of  the  Instantaneous  frequency  Is  obtained  by  taking  the  reciprocal  of  the 
time  between  two  crossings,  then  results  will  be  available  at  an  Irregular 
rate.  On  the  other  hand.  If  we  count  N  crossings  during  an  Interval  of  time  T, 
thereby  maintaining  a  constant  sample  rate,  the  frequency  estimate,  N/T,  is  an 
average  of  instantaneous  frequency  during  time  T.  In  addition,  the  error  of 
the  estimate  will  be  plus  or  minus  one  count,  or  *1/T. 

In  the  present  readout  crossings  are  counted  to  Insure  a  uniform  sampling 
rate  and  the  Interval  T  Is  kept  as  small  (.1  second)  as  possible  In  order  to 
minimize  the  amount  of  averaging,  while  still  maintaining  compatlbllty  with 
the  maximum  data  rate  capability  of  the  acquisition  computer.  The  *1/T  error 
problem  Is  handled  with  the  use  of  a  phase  lock  loop  (PLL)  circuit  (Gardner, 
1979)  that  multiplies  the  difference  frequency  by  a  factor  of  3200,  reducing 
the  scale  factor  of  the  acceleration  and  gravity  Information  to  (15000/3200) 
or  about  4.8  mllllgal/Hz.  The  output  of  the  PLL  is  counted  It  a  specially 
designed  frequency  counter  Incorporated  In  the  readout;  the  resulting  counts 
are  made  available  to  the  acquisition  system.  In  this  way,  the  1/T  (10  Hz  at 
Ta.lsec)  error,  with  the  reduced  scale  factor.  Is  decreased  from  ( 15000x( 1/T) ) 
to  (4.8x( 1/T> ) ,  or  about  48  mllllgals/data  point. 

Clearly,  this  48  mill Igal  error  (at  a  10  Hz  sample  rate)  Is  In  Itself 
unacceptable.  The  frequency  counter  In  the  readout,  however.  Is  designed  so 
that  no  counts  are  missed  between  adjacent  counting  Intervals.  This  enables 
the  data  to  be  processed  such  that  measurement  error  Is  reduced  In  a  way 
unique  to  systems  In  which  Information  Is  encoded  using  frequency.  Unlike 
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conventional  acquisition  systems  (where  Information  Is  available  In  a  sampled 
voltage,  and  errors  In  measurement  are  statistically  Independent),  the  sum  of 
M  data  points  in  the  present  system  will  still  retain  an  ambiguity  of  only 
plus  or  minus  one  count  so  that  the  effective  error  is  reduced  to  1/MT. 
Sampling  at  10  Hz  enables  the  system  to  measure  vehicle  accelerations  to  a 
resolution  of  48  mllllgals  with  a  minimum  of  averaging,  while  still  allowing 
measurement  to  an  accuracy  of  48/M  mllllgals  with  Increased  averaging.  The 
averaging  Implemented  by  summing  counts  Is  a  low  pass  filter  with  an 
approximate  cutoff  at  1/MT  Hz.  In  the  processing  of  data  from  the  test  cruise, 
a  value  of  M=100  was  used  so  that  the  effective  filter  cutoff  was  l/(100x.l) 
or  .1  Hz.  with  an  error  of  48/100  or. 48  mllllgal.  Since  the  signal  of 
Interest,  vertical  gravity,  will  appear  In  frequency  bands  far  below  .1  Hz  In 
both  sea  and  air  vehicles,  the  averaging  does  not  affect  required  performance. 

Along  with  multiplication  of  the  difference  frequency,  the  PLL  can  also 
provide  analog  filtering  of  the  Information  encoded  In  the  FM  modulation.  In 
previous  seagoing  systems  (Bowln  et  al ,  1972),  the  main  bulk  of  low  pass 
filtering  necessary  for  obtaining  gravity  from  the  VSA  data  was  done  In  the 
PLL,  with  a  cutoff  frequency  of  .003  Hz.  The  raw  gravity  signal  was  then 
available  In  real  time.  In  the  present  readout,  since  the  system  Is  designed 
with  airborne  gravity  measurement  In  mind,  a  minimum  of  filtering  Is  done  In 
the  PLL,  so  that  more  precise  digital  filtering  techniques  can  be  applied  to 
the  data  and  so  that  more  Information  Is  available  concerning  vehicle 
accelerations.  Since  airplanes  will  travel  at  higher  speeds,  the  frequency 
band  containing  desired  gravity  signals  will  extend  above  .003  Hz  and  will 
overlap  bands  In  which  vehicle  accelerations  are  prevalent.  Using  Kalman 
filtering  techniques,  a  precise  ancillary  navigation  time  series  together  with 
the  VSA  data  can  separate  these  accelerations  from  the  gravity  signal,  and 
Eotvos  corrections  can  also  be  performed.  By  setting  the  PLL  analog  filter 
cutoff  at  1  Hz,  so  that  negligible  frequency  components  appear  beyond  5  Hz  (In 
keeping  with  the  Nyqulst  criterion  for  a  10  Hz.  sampling  rate)  all  frequency 
components  of  both  gravity  and  vehicle  accelerations  below  1  Hz  are  passed, 
and  the  present  system  should  be  well  designed  for  use  In  aircraft. 

In  the  test  cruise,  unfortunately,  a  major  component  of  the  PLL  had  Inade¬ 
quate  dynamic  range  for  maintaining  phase  lock  for  accelerations  above  .05G, 
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while  the  actual  sea  state  accelerations  In  the  passband  of  the  PLL  filter  were 
on  the  order  of  .1G.  A  simple  component  replacement  was  all  that  was  needed  to 
solve  the  problem,  but  the  part  was  unavailable  on  the  cruise.  One  of  the  data 
tapes,  however,  was  obtained  during  a  period  of  relatively  calm  seas;  a  report 
follows  on  the  processing  techniques  applied  to  this  data  set. 

Processing  Sequence 

Figure  5  Is  an  Illustration  of  about  9  minutes  of  raw  data  obtained  from 
the  readout.  As  mentioned  above,  the  data  was  sampled  at  10  Hz.  with  an  Initial 
resolution  of  48  mil  11  gal  s/count,  so  that  the  overall  vertical  span  of  the 
chart  represents  about  .1  G  or  100000  mil 1 1 gal s .  The  first  minute  or  so  of  data 
Illustrates  the  characteristics  of  the  data  when  the  loop  was  out  of  lock,  and 
the  time  when  lock  occurred  Is  denoted.  This  point  occurred  Immediately  after  a 
course  change,  so  that  the  vessel  was  riding  more  smoothly  and  lock  was  obtain¬ 
able.  In  Figure  6,  about  80  seconds  of  a  horizontally  expanded  section  of 
locked  data  Is  shown  which  details  the  accelerations  due  to  sea  state  as 
measured  by  the  VSA.  Again,  the  total  vertical  scale  encompasses  a  range  of 
about  .1G,  the  resolution  If  48  mllllgal,  and  we  observe  that  the  average 
period  of  wave  motion  was  about  5  seconds  (.2  Hz). 

The  first  step  In  processing  this  data  Involved  the  use  of  a  "moving 
average"  operation  which  simply  replaced  each  point  In  the  data  by  the  sum  of 
100  points  about  It.  In  tills  way,  as  we  have  explained,  resolution  was 
Increased  to  about  .48  mllllgal s/count,  and  an  effective  filtering  of  the  data 
at  a  cutoff  of  about  .1  Hz  was  a  byproduct.  This  "filter"  can  be  viewed  as  a 
finite  Impulse  response  (FIR)  filter  with  100  coefficients,  all  of  which  are 
unity.  At  this  point,  the  output  rate  of  this  filter  was  still  at  10  Hz. 

The  actual  Implementation  of  this  operation  was  done  In  the  frequency 
domain,  to  obtain  Increased  computational  efficiency.  Data  was  taken  In  3000 
point  segments,  an  FFT  was  applied,  the  product  of  the  result  and  the  FFT  of 
the  FIR  coefficients  was  formed,  and  an  Inverse  transform  to  this  product  was 
computed  and  combined  with  results  from  other  segments  using  the  "overlap/add 
method"  (Oppenhelm  and  Schafer,  1975). 

Since  these  data  were  obtained  at  sea  the  gravity  signal  was  In  a  very  low 
frequency  band  and  corrections  for  vertical  motions  of  the  vessel  were 


unnecessary  (they  average  to  zero  over  long  time  periods).  The  large  volume  of 
data  (about  165000  data  points  for  4.5  hours  for  a  distance  covered  of  about 
100  km)  was  unnecessary,  and  a  decimation  of  the  data  by  a  factor  of  100  (l.e. 
a  new  data  rate  of  10  sec)  was  desired.  Before  this  was  done,  however,  a 
second  low  pass  filtering  operation  was  performed  so  that  the  resulting 
decimated  signal  would  satisfy  Nyqulst  criteria  (Crochlere  and  Rablner,  1981). 
A  FIR  filter  was  desired  since  linear  phase  characteristics  would  Insure 
preservation  of  waveforms  such  as  the  gravity  signal  that  are  within  the 
passband.  305  filter  coefficients  were  calculated  with  an  algorithm  that 
maximizes  flatness  In  the  passband  (Parks  and  McClellan,  1972)  with  a  -35  db 
stopband  beyond  .05  Hz.  This  filter  was  again  Implemented  with  the  overlap/add 
technique.  After  selecting  one  out  of  every  hundred  points  of  the  output  of 
this  filter,  the  data  set  was  In  a  more  manageable  size,  and  a  final  50 
coefficient  filter  with  an  effective  cutoff  at  .005  Hz.  was  used  for  the  final 
raw  gravity  result. 

Figure  7  Illustrates  4.5  hours  of  processed  raw  gravity  data  obtained  with 
these  techniques.  The  vertical  scale  now  encompasses  a  total  span  of  400 
mllllgals  or  .00041G.  There  are  two  time  periods  In  which  data  Is  obviously 
Invalid.  Figure  8  Is  a  profile  of  the  east  component  of  ship  velocity  during 
this  time  period,  and  we  see  that  the  patches  of  Invalid  data  occur 
Immediately  before  and  after  significant  changes  In  the  vessel's  course.  In 
the  quieter  sections  of  Figure  7,  abrupt  changes  In  level  are  seen  to  occur  at 
times  which  also  correspond  to  changes  In  east  velocity,  an  Indication  of  the 
presence  of  the  Eotvos  effect  on  gravity  measurements.  An  Eotvos  correction 
was  applied  to  the  raw  gravity,  using  the  east  velocity,  and  the  platform 
latitude  data  available  on  the  initial  data  tapes.  The  results  are  shown  In 
Fig.  9,  In  which  we  observe  that  the  step  changes  have  been  removed  by  the 
corrections.  The  good  quality  of  the  results  of  the  Eotvos  correction,  which 
were  based  on  navigation  data  from  an  Independent  source,  confirms  that  the 
processed  data  Is  Indeed  a  profile  or  free-air  gravity. 

We  have  shown  that  the  raw  YSA  data  obtained  with  the  new  platform  and  PLL 
characteristics,  is  capable  of  measuring  sea  state  accelerations  to  an 
accuracy  of  48  mllllgals  (.00005G).  This  capability  Is  required  for  work 
aboard  aircraft  In  tandem  with  a  second  set  of  acceleration  data  obtained  from 
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Figure  9  Filtered  acceleration  data  summed  with  Eotvos  correction  estimated  from  velocity 

east.  Same  scales  as  for  figure  7.  Note  that  the  discontinuities  at  course  changes 
in  figure  7  are  not  evident  here. 


a  precise  navigation  stream,  so  that  gravity  data  could  be  sorted  out  in 
certain  frequency  bands.  At  the  same  time,  after  suitable  processing, 
accelerations  in  frequency  bands  near  DC  (such  as  longer  wavelength  gravity 
components)  may  be  measured  to  an  accuracy  of  at  least  .48  milligal. 

A  second  test  cruise  is  planned  for  late  January,  1983.  In  order  to 
obtain  a  complete  accurate  gravity  profile  for  the  entire  length  of  the 
cruise,  based  on  the  results  mentioned  above,  we  are  concentrating  on  the 
following  problems: 

1)  Adjust  frequency  response  in  servo  amplifiers 

for  better  performance 

2)  Improve  the  platform  temperature  regulation 

3)  Add  an  error  handling  routine  to  platform  computer 

4)  Calibrate  the  gyros  and  accelerometers 

5)  Improve  the  "navigate  mode"  algorithm 

6)  Redesign  the  phase-lock  loop 
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